Possible direct method to determine THE radius of a star 
from the spectrum of gravitational wave signals 
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We computed the spectrum of gravitational waves from a dust disk star of radius R inspiraling 
into a Kerr black hole of mass M and specific angular momentum a. We found that when R is 
much larger than the wave length of the quasinormal mode, the spectrum has several peaks and 
the separation of peaks Aw is proportional to R~^ irrespective of M and a. This suggests that the 
radius of the star in coalescing binary black hole - star systems may be determined directly from 
the observed spectrum of gravitational wave. This also suggests that the spectrum of the radiation 
may give us important information in gravitational wave astronomy as in optical astronomy. 
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A number of worldwide projects on the direct detection 
of gravitational waves using laser interferometers such as 
LIGO, VIRGO, GEO600, TAMA300, and LISA as well 
as several projects using resonant type detectors are in 
progress . One of the most important sources of grav- 
itational waves for these detectors is coalescing binary 
black hole (BH) - star systems. In the inspiral phase, 
the radius of the star has little effect on gravitational 
waves so that the star can be treated as a point particle. 
Comparing the waveform of the inspiral phase with a the- 
oretical template using matched filtering techniques, we 
may determine each mass and spin of the star and BH, 
respectively 0. While, in the final merging phase, the 
effect of the radius is important and some information 
on the radius of the star may be extracted from gravita- 
tional wave signals. However, this merger phase is much 
less well understood than the inspiral one. 

In this Letter, as a model of this final merging phase, 
we computed the spectrum of gravitational waves from 
a dust disk star of radius R inspiraling into a Kerr BH 
of mass M and specific angular momentum a. We found 
that when R is much larger than the wave length of the 
quasinormal mode (QNM) of Kerr BH, the spectrum has 
several peaks and the separation of peaks Aw is propor- 
tional to R^^ irrespective of M and a. This suggests 
that the radius of the star in coalescing binary BH - star 
systems may be determined directly from the observed 
spectrum of gravitational waves without any assumption 
to equation of state. 

First, let us consider gravitational waves from a par- 
ticle plunging into a Kerr BH of mass M and specific 



*Vallisneri [g] recently argued that the radius of the star 
may be determined if the frequency of gravitational waves at 
the tidal disruption is obtained observationally. However, his 
method needs some information about the equation of state 
P to know the frequency at the tidal disruption. 



angular momentum a. The Fourier component of the ra- 
dial wave function of gravitational waves Ximu{f*) obeys 
the Sasaki-Nakamura equation (Eq. (2.28) of Q) given 
by 
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where r*, Simuj{r*), F{r*) and U{r*) are the tortoise co- 
ordinate of Kerr BH, the source term from the test par- 
ticle of mass ^ (Eq. (2.29) of Q) and two potential func- 
tions (Eqs. (2.12a) and (2.12b) of [§), respectively To 
solve Ximujir*) using a Green function method, we need 
two independent homogeneous solutions whose boundary 
conditions are given by 
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where k ^ uj - ma/[2{M ^/AP - a^)]. Then, the inho- 
mogeneous solution to Eq. ([|) becomes 
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where W is the Wronskian, 
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The asymptotic behavior of the radial wave function 
Xi„iu>ir*) is given by 
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Then, the energy spectrum and waveform of gravitational 
waves (Eq. (3.6) and (3.10) of ||]) are given by 
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where r is the coordinate radius from the center of BH, 
-'i^hn spin - 2 weighted spherical harmonics, cg is 

a constant given in 

Next, we consider gravitational waves from a dust disk 
We note that a geodesic in the equatorial plane 
in Kerr BH is characterized by the specific energy {E) 
and angular momentum (L^). Let t = T{r) and (j) — 
$(r) express an orbit of the geodesic for given E and 
Lz- Then another geodesic for the same E and Lz is 
expressed hy t = T(r) + ct and 4> = ^{r) + c^ where Ct 
and C0 are constants. Specifically the geodesic, whose 
location aX t — T{rQ) is r = r.; and = (/>i, is expressed 
by t = T{r) + T(ro) - r(r,) and = $(r) + 0, - $(r,). 
It is possible to set a number of particles to form a disk 
of radius R whose center is r — and (j) ~ *i'(''o) at 
t = T(ro). 

Each particle in the disk emits the same gravitational 
waves as the particle in the center of the disk except that 
the phase in t and (j) is different from the central particle 

. Then the amplitude of the gravitational waves from 
the dust disk star is computed as 
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where fmui and A^i^^^^'^'^^^ are a form factor and the am- 
plitude of gravitational waves at infinity for a single 

particle case, respectively. S is the normalization fac- 
tor given by 

rra+R 

S^2 drr<j)o{r). (9) 
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In this Letter, we only show the result for vq = IQM, 
a/M — 0.9, E — 1 and Lz ~ 2M. In this case, the orbit 
rotates about vr angle around the BH within the radius 
of r = rg, and the velocity a± r — 3Af, the location 
where the particle might radiate the characteristic QNM 
frequency, is u ~ 0.48c, where c is the speed of light. As 
for the choice of the parameter tq = lOAf , we can refer 
the result of the relativistic Roche limit (Eq. (5.1) of |^) 
as 
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FIG. 1: Energy spectrum of gravitational waves from a dust 
disk star moving on an equatorial plane in Kerr spacetime 
whose radius is set up at ro = IQM for the case of a/M = 0.9, 
Lz/M = 2 ((a) R/M = (test particle), (b) R/M = 1.56, 
(c) R/M = 5.88). We only show I = 2 mode. Solid, dashed, 
dash-dotted, dotted, and dash-three dotted lines denote the 
case of m = —2, —1, 0, 1, 2, respectively. 



M < 4.64Af, 











1/2 



V 10km J 



3/2 



6M\ 



(10) 



where Ms and Tc are the mass of the star and the Roche 
radius, respectively. For appropriate choice of parame- 
ters rp — lOA'f corresponds to the Roche limit. 

Fig. (a) and Fig. || (a) show the characteristics of 
gravitational waves from a test particle plunging into a 
Kerr BH of ajM = 0.9. Fig. |l| (a) shows the spectrum 
for ^ = 2 and Fig. |2| (a) shows the waveform observed 
from the equatorial plane at infinity. From Fig. ^ (a), 
we see that most of gravitational waves are radiated near 
the frequency of the QNM since the QNM frequency for 
a/M = 0.9 is Re (A/wqnm) = 0.7 for Z = m = 2, while 
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the imaginary part of QNM (Im (Mwqnm)) in this case 
is 0.07 (Fig. la of @ and Fig. 3 (c) of §). Fig. | (a) 
shows clearly that there exists a typical ringing tail at 
late time (Fig. 4 of @ and Fig. 4 of (l|]). Fig. | (b), 
(c) and Fig. |^ (b), (c) show the spectra and the wave- 
form from the circular dust disk star inspiraling into a 
Kerr BH with a/M = 0.9. The radius of the disk is set 
up as R/M = 1.56 for Fig. |l| (b) and Fig. | (b) while 
R/M = 5.88 for Fig. (c) and Fig. | (c). For the 
case of Fig. |^ (b), the peak of the spectrum is located 
at Mlo ^ 0.5 which is quite different from the QNM fre- 
quency Re (Mwqnm) ^ 0.7 although the ringing tail can 
be seen in Fig. |2| (b). When we turn to look at Fig. 
|l| (c), there are several peaks in the spectrum and the 
ringing tail is very low or absent in the waveform of Fig. 
§ (c). This is a new finding. The reason for this behav- 
ior is as follows. The energy spectrum {dE / duj)'^^^'' of 
gravitational waves from the disk is expressed as 
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where (d£'/(ia;)|^^'"^'°'' is the spectrum from the single 
test particle. The spectrum from the test particle has 
only one peak at the frequency cjqnm (Fig. [| (a)) so 
that the square of the form factor | fmui P is responsible 
for this behavior. To confirm this, we show | fmu P for 
R/M = 1.56 (Fig. I (a)) and R/M = 5.88 (Fig. | (b)) 
and find that | fmuj P is responsible for the behavior of 
the spectrum in Fig. |l| (b) and (c). The existence of 
several peaks in Fig. || (a) and (b) can be understood by 
the approximate estimation of fmu as 
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which agrees quite well with the numerical results of Fig. 
§ (a) and (b). Equation (p^ ) suggests that the separation 
of peaks of the spectrum Alo may be in proportion to 
R^^. In Table |, we show Alu for various R/M and we 
found that 
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where a constant C is close to 1 in the present case 
the physical unit, Av = Auj/(2n) is given by 
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assuming that there are 4 peaks in the region of G 
[0, wqnm] in the spectrum. Therefore for neutron stars 
and white dwarfs, the frequency band is within some laser 
interferometers and resonant type detectors. 
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FIG. 2: Waveform of gravitational waves from a dust disk 
star moving on an equatorial plane in Kerr spacetime whose 
radius is set up at ro = lOM for the case of a/M — 0.9, 
L^/M = 2 ((a) R/M = (test particle), (b) R/A4 = 1.56, 
(c) R/M — 5.88). We set the angle of the observer at infinity 
as 9 = 7r/2, (p = 0. We only show I — m = 2 mode to exclude 
m mode coupling effect. 



Equation (|Ij) strongly suggests that from the ob- 
served value of Aoj in the spectrum of the gravitational 
wave signals, we may determine the radius of the star if 
R/M ^ 1.6. From the physical point of view, the form 
factor fmuj expresses the phase cancellation effect of grav- 
itational waves 1^, so that if i? < I/wqnm gravitational 
waves from the dust disk star is essentially regarded as 
that from the test particle. [| This suggests that our 



^ We also confirmed that this phase cancellation can also 
be seen in the 3D dust star assuming that each test parti- 
cle moves the constant polar angle with no orbital angular 
momentum. 
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proposal is valid only if the inspiraling star is tidally dis- 
rupted by the BH finally. It is instructive to note here 
that in reality the form factor may bring us the informa- 
tion about the form of the source, radius. 

Here, we only showed the results for a special set of pa- 
rameter such as ro = lOM, a/M = 0.9, E = 1, = 2M. 
To check the robustness of our proposal, we also calcu- 
lated the spectra from the dust disk star inspiraling into 
Kerr BH for the wide range of parameters {Lz : spiraling 
case, < a/M < 0.9) and found that there exist some 
peaks in the spectrum for R ^ I/wqnm- In these cases 
we confirmed the relation R 1/ Aw. In this Letter, we 
did not take into account of the pressure and the self- 
gravity of the star so that it is urgent to confirm our 
proposal by full 3D numerical simulations including the 
determination of the value of the constant C in Eq. (14). 
If the relation of Eq. (|l4|) is confirmed irrespective of the 
equation of state, Eq. (|14|) can be adopted as one of the 
direct methods to determine the radius of the star from 
gravitational waves. In any case, it is quite possible that 
the spectrum of the gravitational waves may give us im- 
portant information in gravitational wave astronomy as 
in optical astronomy. 



TABLE I. Comparison with the characteristic length from 
the energy spectrum of gravitational waves to the radius of 
the disk. R denotes the coordinate radius. 




FIG. 3: Form factor of the dust disk star moving on an 
equatorial plane in Kerr spacetime whose radius is set up at 
ro = lOM for the case of a/M = 0.9, L,/M = 2 ((a) R/M = 
1.56, (b) R/M = 5.88). We only show I = m = 2 mode, 
because m = 2 is a dominant mode for above parameters in 
1 = 2. It is clear that the form factor is responsible for the 
spectra in Fig. 1. 
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